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^* Abstract. Hydrodynamics of superfluids can be described by formally dividing the fluid into 

a normal fluid and a superfluid part. In color-flavor locked quark matter, at least one superfluid 

component is present due to spontaneous breaking of baryon number conservation, and an additional 

one due to the breaking of strangeness arises once one takes into account kaon condensation. We 

rS^ show how such a two-component description emerges from an underlying scalar field theory which 

can be viewed as an effective theory for kaons. Furthermore, the occurring hydrodynamic quantities 

ft in the low-temperature limit are related to the microscopic parameters provided by the Lagrangian 

which closes the gap between field theory and hydrodynamics, which are important for astrophysical 

calculations. 

> 

1 Introduction 

The study of quantum chromodynamics (QCD) at low temperatures and high den- 
sities is relevant for fundamental as well as applied astrophysical questions. On the 
-^ fundamental side, the phase diagram of QCD at high and intermediate densities might 

t-h be very complicated and its study is highly non-trivial. While heavy-ion collisions and 

£> lattice calculations are powerful tools to probe the phase structure in a regime of high 

temperatures and low densities, their applicability at higher densities is very limited. 
However, future accelerator facilities such as NICA might provide some insight (see 
for example Ref. PQ). 

Reliable information can be obtained in a region of asymptotically high densities, 
where QCD behaves like a free field theory and perturbative calculations are possible. 
The ground state in this region of the phase diagram is a color superconductor where 
quarks of all color and flavor form cooper pairs which was termed color-flavor locking 
(CFL) [2]. Among many interesting features of CFL that could be discussed at this 
point, two are particularly important for the following analysis: CFL spontaneously 
breaks chiral symmetry which leads to an octet of Goldstone bosons and furthermore, 
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CFL spontaneously breaks baryon conservation U(1)b- Strictly speaking, perturba- 
tive calculations of for example the magnitude of the superconducting gap, are only 
valid at (rather exotic) values of the chemical potential of about // ~ 10 8 MeV and it is 
thus questionable whether CFL will persist all the way down to the phase boundary of 
nuclear matter. Going down in density, the increase in the mass of the strange quark 
will act as an external stress on the highly symmetric pairing pattern and systematic 
studies show |3], that CFL will most likely react by developing a kaon condensate. 
The corresponding ground state CFL-K is then subject to two symmetry breaking 
patterns: The breaking of chiral symmetry and baryon conservation in CFL and the 
breaking of strangeness conservation due to the kaon condensate. 

SU(3) L ®SU(3) R ®SU(3) C ®U(1) B -> SU{3) L+R+C ®Z 2 , SU(3) L+R+C 2 U(l) s -> 1 . 

(1) 
One should note at this point that weak interactions explicitly violate strangeness 

conservation and thus U(l)s is not an exact symmetry to begin with (we shall come 
back to this point in the outlook). Apart from CFL-.K' , there are many other possible 
candidate ground states of QCD at intermediate densities [I], but the absence of 
controlled experiments in laboratories or reliable models to perform calculations make 
it very difficult to decide which phases are realized in nature. 

The study of compact stars as the only "laboratory" where such intermediate den- 
sities are realized could prove to be very useful. The maximum value for the chemical 
potential inside a compact star is estimated to be high enough for deconfined quark 
matter to be conceivable (fi s t ar i$ 500 MeV) . Usually, most insights into the interior 
of compact stars can be obtained from the characteristics of cooling and neutrino 
emissivity as well as from mass-radius relations. Other observable phenomena such 
as pulsar glitches or r-mode instabilities require a hydrodynamic description. In order 
to provide a framework for such calculations, it is important to understand how an 
effective hydrodynamic description emerges from the underlying microscopic physics 
discussed above. Moreover, different (equivalent) formulations of superfluid hydrody- 
namics exist in the literature and should be connected to each other. These were the 
central aims of Ref. j5] whose main results shall be reviewed here. 

2 Superfluid hydrodynamics from field theory 

The spontaneous breaking of a U(l) symmetry, be it the breaking of U(1)b or U(l)s 
discussed above, is known to be the fundamental mechanism for superfluidity. We will 
focus here on kaon condensation (and thus on the breaking of U(l)s), but similar con- 
siderations hold for essentially any system with a spontaneously broken symmetry. 
The existence of a condensate together with the absence of elementary excitations 
which could dissipate energy allow for frictionless transport of the associated U(l)s 
charge at low temperatures. This might be puzzling at first since it is also well known 



that the spontaneous breaking of a global symmetry leads to existence of massless 
excitations (Goldstone bosons). At least at zero temperature, this apparent contradic- 
tion can easily be resolved: since the low-energy dispersion relation of the Goldstone 
mode is linear in momentum, it can only be excited beyond a certain critical velocity 
which is given by the slope of the linear part of the dispersion. At nonzero tem- 
perature, the situation is more complicated as thermal excitations of the Goldstone 
mode ("phonons") are present for any superffuid velocity. As a consequence, in the 
translation into hydrodynamic equations, we will have to consider both ingredients, 
condensate and phonons, and they also appear coupled to each other. 

The hydrodynamic framework capable of describing superffuids, which we will 
derive now from microscopic physics, was first set up by Landau [6] mainly for the 
(non-relativistic) purpose of describing superfluid helium. The basic idea is to for- 
mally divide the fluid into a superfluid and normal fluid part. The total mass density 
(later in the relativistic context to be replaced by charge density) is then given by 
P = pn + Ps where at zero temperature only the superfluid density is present and 
above the critical temperature only the normal fluid density. From this point of view, 
the hydrodynamic description of CFL with kaon condensation is quite complex: due 
to the breaking of U(1)b and U(l)s two superfluid components are present such that 
in total we expect to be dealing with four different fluid components. In order to 
derive the two component hydrodynamics for the kaon part, we start by analyzing a 
complex scalar ip 4 model, which can be interpreted as an effective theory for kaons 

0- 

2.1 Zero temperature 

Based on Ref. [5], we start with our derivation at zero temperature, where only the 
superfluid component is present. As mentioned at the end of the last section, our 
starting point will be a ip A theory where the complex doublet field can be interpreted 
as (<^,<^*) = (K ,K ): 

C = <9^dV - ™ 2 M 2 - A M 4 • (2) 

After introducing the condensate, 

cp ->• (<p) + fluctuations , (<p) = ^- e ^ {x) , (3) 

v2 

we simplify calculations by additionally assuming a constant modulus p of the con- 
densate. The equations of motion then imply [3] that also d^ip is constant and the 
tree-level potential turns into: 

u = -c = _^!z_p!l! > o^^/d^p. (4) 



Field-theoretic definitions for the Noether current and stress energy tensor can be 
matched to the corresponding hydrodynamic expressions: 

f = ^ = mr) =a ^^r-- (5) 

9 S(*/=a£) a 2 - m 2 
T"" = (e s + P,)« - gTP. = = W x - ] = WW* > cT£.(6) 

Charge density n s and the flow velocity of the superfluid v% can than easily be ob- 
tained: 

a 2 — m 2 d^ip 



n s = vW" = o- r , < = — - , (7) 

A U 

as can be energy and pressure density of the superfluid (e s and P s ) by making use 
of proper projections on T^ v . It is important to notice that our simplification of 
d^ip being constant translates into a uniform superflow velocity. Finally we include 
thermodynamics into our considerations which allows us to identify the chemical 
potential in the superfluid rest frame, fi s = a. The Lorentz invariant quantity o can 
also be expressed in terms of the chemical potential in the normal-fluid rest frame 
H = doip: with the above definition of the superfluid velocity, a usual Lorentz factor 
arises in this relation, 

a = ^l-v 2 . (8) 

As we can see, both chemical potential and flow velocity of the superfluid evolve solely 
from rotations of the phase of the condensate: rotations of the phase around the U(l) 
circle per unit time create the chemical potential and the rotations per unit length 
gives rise to the superflow velocity. 

2.2 Finite temperature 

In order to introduce finite temperature T, we use the simple one-loop effective action 
with the tree-level potential U as introduced before and the inverse propagator S -1 : 

V 1^ S~\k) , ( -k 2 + 2(a 2 -m 2 ) 2ik • ty \ 

1 ~~T 2^ ~T 2 ~' lj ~V -2ik-dil> -k 2 )' 

(9) 
We work at small coupling A <$; 1 and restrict ourselves to small temperatures (much 
below the critical temperature and T <C /x). Furthermore, any temperature depen- 
dence of the condensate p has been neglected. These approximations allow us to ob- 
tain analytic results (generalizations are discussed in the outlook). Finite-temperature 
calculations are carried out in the Matsubara formalism, i.e. k = (ko, k), ko — —iu) n 
with bosonic Matsubara frequencies u n = 2nnT. For small momenta, the dispersion 



of the Goldstone mode obtained from the poles of the propagator can be written as: 

e l : = c 1 (x)\k\ + ^-\k\ 3 + ..., (10) 

and obviously depend on x = cos?? where i? is the angle in between k and the ho- 
mogeneous background superflow v s . For explicit results for c\ and c 2 please refer to 
Ref. 0. 

At finite temperature, a hydrodynamic interpretation of field-theoretic results is 
much more challenging. We now have to introduce an additional "normal" current 
related to the elementary excitations which automatically complicates thermodynam- 
ics: since our two homogeneous flows imply two rest frames where either the normal 
or the super current is vanishing, it is impossible to define a frame where the pressure 
is isotropic. We rather have to deal with pressure densities normal and longitudinal 
to the non vanishing current (T±,T\\) in the respective rest frame, but which pressure 
is the correct one to appear in the thermodynamic relation? 

2.2.1 Generalized thermodynamics and entrainment 

The most rigorous answer to this question is contained in the canonical formalism 
introduced by Carter [9] which provides an extension of usual thermodynamics by 
making use of a generalized pressure \l/ and a generalized energy density A. The set of 
variables in this framework differs from the one used in the two-fluid formalism of Lan- 
dau: instead of formally dividing the Noether current into normal and supercurrent, 
the canonical framework is based on either the conserved charge and entropy currents 
j M and s M (the latter is conserved only in the non-dissipative case) or their thermody- 
namically conjugate momenta d^ip and 9^. Lorentz covariance requires A and \l/ to be 
functionals of invariants only (i.e. A = A[j 2 , s 2 , j^s 1 *], \& = ^[a 2 ,9 2 ,d /Ji ip9 IJ ']) . Here 
the definition of a from above has been used. The relation of d^ip to the chemical 
potential has already been discussed in the previous section and as we show later, 6> M 
can be related to the temperature. It might seem puzzling at first to generalize scalar 
quantities such as temperature to a four-vector, but the actually measured tempera- 
ture is of course obtained by contraction with the velocity four-vector pointing to the 
rest frame of interest. We can now obtain the correct equations of motion by varying 
either the generalized energy density with respect to the conserved currents or the 
generalized pressure with respect to the conjugated momenta: 

dK = d^df + M ds" , d& = 3^(3^) + s M d0" . (11) 

Switching between these two descriptions is identical to a change from a Hamiltonian 
to a Lagrangian framework, A and \l/ are connected by a Legendre transformation. 
With all these ingredients at hand, we can now write down a generalized thermody- 
namic relation: 

A+qr=j-dip + s-9. (12) 
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Applying the chain rule to Eq. (11) allows us to write: 

<9A <9A 

&>$ = — = Bf + As^ , 0" = — = Af + Cs^ , (13) 

with 

<9A E _ 9 dA r-2 dA 

d(j ■ s) ' " dj 2 ' ' ds 2 

The coefficient A appears in the relations for d^ifi as well as 9^. It is called "entrain- 
ment coefficient" which reflects the fact that any change in either of the two conjugate 
momenta will automatically affect both currents (and the other way around) such that 
the currents are coupled to each other. Finally, the stress-energy tensor is now given 

by 

T^ = -g^V + fd u ^/j + s»6 v . (15) 

Plugging in Eqs. (13), one can see right away that T^ u becomes symmetric in the 
Lorentz indices. 

It remains to relate the canonical formalism to the original two fluid picture of 
Landau (a relativistic generalization of which can for example be found in [TOl [TT]). 
After introducing a normal-fluid current, the straightforward extension of our zero- 
temperature expressions for Noether current and stress-energy tensor are 

T^ = (e n + P n )uy n + (e s + P s )d^d^/a 2 + (P n + P s )g^ , (16) 

f = n n u£ + n s d^/a . (17) 

Here, we have introduced the flow velocity of the normal fluid u£, as well as the 
corresponding energy and pressure densities e n and P n . The superfluid velocity is still 
given by Eq. (7) . The coefficients A, B and C prove to be very useful in translating 
one formalism into the other as all hydrodynamic parameters in either model can be 
expressed in terms of them: 

° A r, ° 2 r, BC-A 2 2 

n s=»' n n = -^s, e s + P s = —, e n + P n = s . (18) 

Note that in (16) and (17) j M and at the same time d^ip appear as parameters such 
that in the strict sense of the canonical framework, this has to be considered a "mixed" 
form. 

2.2.2 Explicit results at finite temperature 

In the terminology of the Landau formalism, it can be shown [5] that all microscopic 
calculations at finite temperature are carried out in the normal fluid rest frame where 
we can identify generalized pressure \1/ and temperature: 

T 
^ = — T , 9q = T in the normal fluid rest frame. (19) 
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Figure 1: Velocities of first and second sound u\ and U2- The homogeneous back- 
ground superflow is aligned along the horizontal axis and the speeds of sound are 
calculated for all angles with respect to the superflow. 



Even though the first relation is in principle frame independent, using our microscopic 
expression (9) automatically restricts us to the normal rest frame. We can now step 
by step derive field-theoretic definitions for all the basic quantities of the two-fluid 
formalism. The coefficients A, B and C are in terms of microscopic variables 



A 



d°ip 



s°j ■ Vip 



V 



B 



(WO 2 



c 



j°d°ipri-j-Vil)s 6 



;CA2 



j-W 



(20) 



with the abbreviation rj = v 2 j d°ip + j ■ Vip. With the help of these, we can easily 
obtain explicit results of for example the normal-fluid density in the low-temperature 
limit, 

4vr 2 T 4 (1 - v 2 ,) 2 48vr 4 T 6 (1 - v 2 ) 2 , „ .. 

(15 + 38^-9^) . (21) 



n,, 



5^3 \i (1 - 3^)3 Tv/3 /i 3 (1 - 3%y 

As it turns out, the normal fluid density is not just the phonon contribution to the 
number density which would yield a T 3 rather than T 4 //z contribution. 

Finally, as an application we present the velocities of first and second sound in the 
background of the superflow. The complete calculation is lengthy and can be found 
in the appendix of Ref. [5]. The final results for T = are shown in Fig. 1. The 
sound velocities depend on the angle between the direction of the sound wave and 
the direction of the superflow. The velocity of second sound decreases significantly 
as the superflow approaches its critical velocity v crit = l/\/3. The impact of small 
temperature corrections on the speeds of sound is also discussed in Ref. [5] . 



3 Outlook 

For the sake of obtaining analytic results, simplifications were made at several points. 
This was a necessary and important first step to gain an understanding of how su- 



perfluid hydrodynamics emerge from the underlying microscopic physics. Restriction 
to a low-temperature regime can be overcome by making use of the two-particle irre- 
ducible effective action pj H2] . All hydrodynamic parameters can then be calculated 
numerically up to the critical temperature (and beyond) and the temperature de- 
pendence of the condensate can be taken into account. Furthermore the explicit 
breaking of U(l)s due to weak interactions can be taken into account by adding a 
small symmetry breaking term to the Lagrangian. 
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